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C 1 ABSTRACT 

a " 

We discuss about the speed of sound squared in the Mass Varying Neutrinos scenario 

co ; 

(MaVaNs). Recently, it was argued that the MaVaNs has a catastrophic instability which 

CN 

J> ' is the emergence of an imaginary speed of sound at the non-relativistic limit of neutrinos. 

o\ : 

As the result of this instability, the neutrino-acceleron fluid cannot act as the dark energy. 
However, it is found that the speed of sound squared in the neutrino-acceleron fluid could be 
positive in our model. We examine the speed of sound in two cases of the scalar potential. 
One is the small fractional power-law potential and another is the logarithmic one. The 
Q ■ power-law potential model with the right-handed neutrinos gives a stable one. 
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I INTRODUCTION 



One of the most challenging questions in both cosmological and particle physics is the nature 
of the dark energy in the Universe. At the present epoch, the energy density of the Universe 
is dominated by a dark energy component, whose negative pressure causes the expansion of 
the Universe to accelerate. In order to clarify the origin of the dark energy, one has tried to 
understand the connection of the dark energy with particle physics. 

In a scenario proposed by Fardon, Nelson and Winer (MaVaNs), relic neutrinos could form 
a negative pressure fluid and cause cosmic acceleration pQ. In this idea, an unknown scalar 
field which is called "acceleron" is introduced and neutrinos interact through a new scalar 
force. The acceleron field sits at the instantaneous minimum of its potential, and the cosmic 
expansion only modulates this minimum through changes in the neutrino density. Therefore 
the neutrino mass is given by the acceleron, in other words, it depends on its number density 
and changes with the evolution of the Universe. The cosmological parameter w and the dark 
energy also evolve with the neutrino mass. Those evolutions depend on a model of the scalar 
potential strongly. Typical examples of the potential have been discussed by Peccei 0. 

The variable neutrino mass was considered at first in [3: , and was discussed for neutrino 
clouds jlj. Ref. [5] considered coupling of the dark energy scalar, such as Quintessence to 
the neutrinos and discuss its impact on the neutrino mass limits from Baryogenesis. The 
MaVaNs scenario leads to interesting phenomenological results. The neutrino oscillations 
may be a probe of the dark energy 00 El- The leptogenesis 0110], the cosmo MSW effect of 
neutrinos ^T] and the solar neutrino [T^nTIj have been studied in the context of this scenario. 
Cosmological discussions of the scenario are also presented [HJ EH HZl UH E3 120] • The 
extension to the supersymmetry have been presented in ref. |2*T| 12*2*] . This scenario is also 
discussed in the context of the texture of the neutrino mass matrix with three families j2Hj- 

Despite many implications of this scenario, ref. [21] showed that this scenario contains a 
catastrophic instability which occurs when neutrinos become non-relativistic. As neutrinos 
become non-relativistic, the speed of sound squared in the neutrino-acceleron fluid turns to 
be negative. As the result of this instability, neutrinos condense into neutrino nuggets, and 
thus cannot act as the dark energy. Here, it is important to say that some non-adiabatic 
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models |TU do not suffer from this instability. However, we have found that the speed of 
sound squared in this fluid could be positive even though neutrinos are not enough relativistic 
and the neutrino-acceleron fluid is adiabatic. In order to realize the positive speed of sound 
squared, a constraint for the scalar potential is required. 

The paper is organized as follows: in Sec. II, we summarize the MaVaNs scenario with 
three families. Sec. Ill presents discussions for the speed of sound in the hydrodynamic 
picture. Sec. IV presents a stable model in the MaVaNs scenario. Sec.V devotes to the 
summary. 

II DARK ENERGY FROM MaVaNs 

In the MaVaNs scenario, one considers a dark energy sector consisting an "acceleron" field, 
4> a and a dark fermion, ip n . This sector couples to the standard model sector only through 
neutrinos. The dark energy is assumed to be the sum of the energy densities of neutrinos 
and a scalar potential for the acceleron: 

Pde = pv + V((f) a ), (1) 

where the potential energy of the acceleron is responsible for the acceleration of the Universe 
and for the dynamical neutrino mass. The energy density for three generations of neutrinos 
and antineutrinos is generally given by 

* = r« 6^, F^^f^^., (2) 

~y T n z Jo e y + 1 

where i denotes three families. 1 

In the scenario, pde is stationary with respect to variation in the neutrino mass. This 
stationary condition is represented by 

dpv dV((p a (m ui )) _ 

3 + 3 - U. {6) 

i=l i=l 

1 In ref. there is a more precise discussion about the Fermi factor in eq.JJJ). For our purpose it 

suffices to use this form. 
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If dJ2 m vi/d(fia 7^ 0, this condition turns to 

T ^f|L + f^). . (4) 
^ dii d(j) a d(p a 

Using the equation of energy conservation in the Robertson- Walker background and the 
above stationary condition, one can get the equation of state parameter w as follows: 

a + l J^SM, (5) 

where 

l^i ?• at. 

h(T) = ^ (6) 

0=1 

The speed of sound squared in the neutrino-acceleron fluid is given by 

c s = — = : (7) 

Pde Pde 

where p is the pressure of the dark energy |2H1- Recently, it was argued that when 
neutrinos are non-relativistic, this speed of sound squared becomes negative in this scenario: 
c 2 s = {d\n.m v / d\n.n v ) < 0, where n v is the number density of neutrinos. The emergence of 
an imaginary speed of sound shows that the MaVaNs scenario with non-relativistic neutrinos 
is unstable, and thus the fluid in the scenario cannot act as the dark energy. However, it is 
found that the speed of sound squared in this fluid can be positive even though neutrinos 
are non-relativistic. Then, a constraint for the scalar potential is required. 

Ill SPEED OF SOUND 

At the non-relativistic limit of neutrinos, the energy density of neutrinos is given by 

3 

p u = ^2m ui n v , (8) 
i=i 

then the stationary condition eq.Q is rewritten as follows: 

= (9) 



Now the dark energy density is given as 

3 

Pde = E m vi n v + V(<p a ). (10) 

i=l 

Using these relations, the speed of sound squared in the neutrino-acceleron fluid is negative: 
c 2 = (d\n.m v / d\n.n v ) < as shown in ref. P2j. In order to study c 2 quantitatively, we start 
with discussing the energy density of neutrinos in eq.(J2J). Taking account that is much 
larger than 1 for non-relativistic neutrinos, the function -F(£«) is expanded in terms of 
as: 

where 

T 3 too (j m P coo dyy± 

h =— ayy a = Jo eV + 1 -6 47 (12) 

v ix 2 Jo ey + 1 ' 2 / °° ^ ' 

Since the first term of the right hand side in eq. (jl(J|) is derived from the first term of the 

right hand side in eq.(JTTJ), the effect of the second term in eq.(JTTJ) should be added to the 

dark energy density as a correction. Due to this correction, it could be that the negative 

speed of sound squared turns to be positive under a condition, which is discussed in detail 

later. 

The dark energy density including the correction term is given as follows: 



Pde = E rn vi h p I 1 + ) 

i=l V Si / 



+ V{<f> a ). (13) 



Then the stationary condition eq.(0]) is described as 



fr[ V & J d(p a d(f) a 
The equation of state is generally given by eq.©. At the non-relativistic limit, it is easy to 
see h(T) = 1. Because of the correction term, h(T) is deviated from 1 as follows: 



Thus, the equation of state is given by 



w + l= 1=1 \ ^-L (16) 

where we omitted the term of (9(l/£ 3 ). Using eqs. (|13|) and 1)1 6j) and the stationary condition 
eq.(j2J), the speed of sound squared in the neutrino-acceleron fluid is described finally as 
follows (see Appendix): 



c 



PDE ~ 

az 

A dm ui A 3 / 5T T dm ^ 

f=i dz 3 fri V Ci H dz 



+ * =i v ^ — (17) 

dn u 3 dr^ 

f=i dz ti dz 

where z is the redshift parameter, z = (T/Tq) — 1 and "0" represents a value at the present 
epoch. We have taken the form differentiated by the redshift parameter instead of the time. 
The first term in the right hand side of eq. (|17|) is the leading term at the non-relativistic 
limit and is negative definite because of dm^/dz < and dn u /dz > 0. The numerator of the 
second term, which is the correction term, is the positive definite. Thus, it is possible that 
the speed of sound squared is positive, in other words, the neutrino-acceleron fluid could be 
stable due to this term. From the eq.([17|h it is easy to see that if the following relation is 
satisfied, the speed of sound squared becomes positive: 

h~{ ^) — 3 — ( ' 

The sign of (? s is determined by the magnitude of dm ui /dz. Since the magnitude of dm ui /dz 
is model dependent, we will examine two cases of the scalar potential for the acceleron, one 
is the power-law potential and another is the logarithmic one in the next section. 



IV MODELS 



In the MaVaNs scenario, one needs a flat scalar potential P]. Therefore, we will discuss the 
speed of sound in cases of the small fractional power-law potential and logarithmic one in 
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this section. The magnitude of dm V i/dz does not only depend on the scalar potential but 
also the coupling between neutrinos and the acceleron. Therefore, we consider two cases in 
eqs. ()2(J|) and (|28|) for this coupling. 

A. The power-law potential 

We take the small fractional power-law potential of the form 

Vfa)=A(^) , (19) 



where parameters A and k are fixed by the magnitude of the dark energy and the stationary 
condition at the present epoch, respectively. 

Three left-handed neutrinos and a sterile neutrino 

We take a Lagrangian of the form 

C = v LoL m a D ^ n + \4> a ipn4>n + h.c. , (20) 

where Ul and ipn are the left-handed and a sterile neutrino, respectively. Since we consider 
three families of active neutrinos, the mass matrix mo is 3 X 1 matrix. Thus the neutrino 
mass matrix is given by 

M. = (21) 

where we assume mf) <C A0 a . After diagonarizing this matrix, we can find mass eigenvalues 
of neutrinos as follows: 

M 2 

m vi = -f-, (22) 

which gives 

M 2 

Using eqs.((H), (fTT?|) and (J2BJ), we have the redshift dependence of the acceleron from the 
stationary condition eq.(jU) as 
-Ak + G(z) 



-, G(z) = 

2\an u T 2 Y.{l/Mf) 
i=i 



\ 



AW + AanlT 2 j^(l/Mj) E M 2 , (24) 

j=l k=l 



where the leading term is taken after expanding dV/d(p a by k. Thus, the redshift dependence 
of neutrino masses is given by 



3 dm vi 
^ dz 



E 



M? 



1 Ml dz 



Aa 2 hlT A 



i=i 



Em! 



i=i 



[z + 1)[-Ak + G(z)]< 



8h u E M i 



i=i 



5[-Ak + G(z)] 

~GTz) ^ ~ 

2ah u T^^/Mf) 
i=i 



(25) 



where the right hand side is negative and its absolute value is the increasing function of k. 
At the present epoch, values of some parameters are given by 



T ~ 1.69 x 10" 4 (eV), A ~ 2.99 x l(T u (eV 4 ), n 



and we take the following typical masses at the present epoch: 



12 x l(T 13 (eV 3 ), 



(26) 



m° vl = 0.0045(eV), m° 2 = O.Ol(eV), m° 3 = 0.05(eV), 



(27) 



which lead to Am 2 n 



atm 



2.4 x 10" 3 (eV 2 ) and Am| un = 8.0 x 10~ 5 (eV 2 ). We numerically 
evaluate the relation (|18j) by putting values of (|26|) and ()27|). We find that k has to be smaller 
than 5.5 x 10~ 5 to satisfy the relation (J18)) . However this value of k is unfavored in the 
phenomenology of the neutrino experiments. The value of k is related with neutrino masses 
through the stationary condition eq. (jl4j) . Actually, the k which is smaller than 5.5 x 10~ 5 
leads to ~ 0(l(T 4 )(eV). In the case of 10" 4 < k < KT 2 , E"i„ ~ 0(l(r 3 ~ KHXeV) 

is expected. Therefore, this model including the relation (|2*2"j) and the scalar potential (|T9"|) 
is unfavored. 

Three left- and right-handed neutrinos and a sterile neutrino 



We add usual right-handed neutrinos to the Lagrangian (|2T)|) 



£ = v Lo m a D ^ n + \(p a ip n ip n + v La M^v Rti + z^Mjf C v v m + h.c. 



(28) 



where z/r is the right-handed neutrinos and both Md and M# are 3x3 matrix. Then, the 
neutrino mass matrix is given as the 7x7 matrix 



M 



I m D M D 

m T D \<p a 
V Ml M R 



(29) 



in the (ul, ip n , vr) basis. We take the right-handed Majorana mass scale to be much higher 
than the Dirac neutrino mass scale, and assume <C \<f) a . Then, the effective neutrino 
mass matrix is approximately given by 

M v = M D M^M T D + (30) 

The first term in the right hand side of eq. (|3U|) is the time- independent neutrino seesaw 
mass matrix, which is denoted by M u , and so it depends on the flavor model of neutrinos. 
In this case, since the neutrino mass matrix reduces to eq. (j21|) at the limit of M u —>■ 0, we 
assume that M v dominates the effective neutrino mass matrix M v . Then we can describe 
generally mass eigenvalues in the first order perturbation as follows [23J: 



M 2 

m vi = m ui + q— f-, (31) 

where c« is a coefficient of order 1 depending on the model of families. 2 Using this relation 
and eq. (jl9|) . we have the redshift dependence of the acceleron from the stationary condition 
eq.jTU) as 



3 



7? 



U \ m li 



XkA • < ffl > 

where the first term of the left hand side in eq. (fT^j) was expanded by the second term of the 

right hand side in eq. (j31|) . Thus, the redshift dependence of neutrino masses is given by 



3 3 / ^nT 2 \ 

kA^M^Mf 3-9- 

i=l j=l \ vj / 



3 dm^ = _ ^ , 

^ r 3 / n rp2^ "I 2 



i=l 9Z 



h u (z + 1) 



aT 2 



2 

k=l \ m vk , 



The magnitude of the first term of the left hand side in eq.([18|) is nearly equal to kA/n u 

at the present epoch. In order to realize the positive speed of sound squared, k has to be 

smaller than 5.19 x 10~ 6 . When we assume that the magnitude of the second term of the 

right hand side in eq.(|31|) is 0.1 percent of the first term, we obtain k = 1.90 x 10~ 6 which 

reproduce observed values of neutrino masses. Therefore, the model including the scalar 

potential (JT§j) and the relation between neutrino masses and the acceleron (}3T|) is favored in 

the MaVaNs scenario and can act as the dark energy. 

2 The values of q are given in a specific flavor symmetry as in ref. \2'.'>\ . 



B. The logarithmic potential 



We take the logarithmic scalar potential of the form 



Bin ^ 



H ) 

-14/ T T 1\ 



B ~ 5.68 x 10 _i4 (eV 



(34) 
(35) 



where we use values of (f2T)|) and (j2ZJ) to fix the value of B. 

We will consider two cases of the coupling between neutrinos and the acceleron as well 
as the case of the small fractional potential. 

Three left-handed neutrinos and a sterile neutrino 

We take Lagrangian in eq.(|2U|). and thus, neutrino masses in eq. ()22|) . The stationary condi- 
tion of eq. (|14j) gives the redshift dependence of the acceleron: 



-B + H(z) 



i=i 



H(z) = 



\ 



B 2 + Aan 2 v T 2 J20-/ M j) E M i 

5=1 



hi 



(36) 



k=l 



and thus the redshift dependence of neutrino masses is given by 



4aVT 4 



^ dz 

i=i uz 



i=i 



E^i 2 



i=i 



(z + l)[-B + H(z)}< 



8iVM, 2 

ti 5{-B + H(z)} 



H(z) 



2an u T 2 J2(l/M? 
i=i 



(37) 



Using values of eqs.(jH)j), (|2*7)l and (|33j). the first term of the left hand side in eq.fTH]) is 

— O(10 _1 ) (eV), however the second term is O(10 _4 )(eV) at the present epoch. Thus the 
speed of sound squared becomes negative, in other words, one cannot build a stable MaVaNs 
model including the form of the neutrino mass like as eq.| 



Three left- and right-handed neutrinos and a sterile neutrino 



We take Lagrangian in eq. (|28j) . and thus, neutrino masses in eq. (J3T| . The redshift depen- 
dence of the acceleron is given as: 



^ 2 / aT r 



i=i 



m li, 



XB 



(3* 



10 



and thus, the redshift dependence of neutrino masses is given by 




5aT 2 



mi 



) 

-i 2 ' 



i=l j=l 



(39) 



hu{z + i) £m, 2 i 



Lfc=i 




aT 2 



) 



Using values of eqs.(J2SJ) and ()27jl. we can estimate the redshift dependence of neutrino 
masses. Since fh V i dominate the neutrino mass, we assume fh u i ~ m^. In this case, the 
first term of the left hand side in eq. (jl8|) is — O(10~ 3 )(eV), however the second term is 
C(10 _4 )(eV). Thus, the speed of sound squared becomes negative as well as the previous 
case. We conclude that the speed of sound squared in the neutrino-acceleron fluid becomes 
negative when neutrinos are non-relativistic for the logarithmic scalar potential in eq. (|34|) 
and neutrino masses like as eqs. ()22j) and (j31j) . 



We have discussed about the speed of sound squared in the neutrino-acceleron fluid, and 
tried to find a stable MaVaNs model in which the fluid is adiabatic. In order to examine c 2 
quantitatively, we have taken two types of the scalar potential for the acceleron. One is the 
small fractional power-law potential and another is the logarithmic one. Furthermore, we 
have studied about two types of models which have different couplings between neutrinos 
and the accerelon. 

In our analysis, models including the logarithmic scalar potential are unstable and cannot 
act as the dark energy because the speed of sound squared in the neutrino-acceleron fluid 
becomes negative. However, models with the small fractional power-law potential can avoid 
this instability. The model including only three left-handed neutrinos and a sterile neutrino 
avoids this instability but does not reproduce the observed neutrino masses. On the other 
hand, the model including the right-handed neutrinos reproduces the observed neutrino 
masses and realizes the positive speed of sound squared. Neutrino masses in this model 
have the time-independent component from the seesaw mechanism, which was assumed to 
be dominant in the effective neutrino mass. Therefore, it is easy to reconcile these neutrino 
masses with observed ones. Due to this time-independent mass, this model becomes viable. 



V SUMMARY 
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Appendix 



The derivation of the speed of sound squared in eq.(JT7J) is presented in this Appendix. 
The energy density for three generations of neutrinos and antineutrinos is generally given 

by 



m, 



i=l 



/ " ./„ <'■< ■ 1 ' 

As neutrinos become non-relativistic, & is much larger than 1. Therefore, the function F(£i) 
is expanded in terms of ^ as: 



dyy 2 



TP 

X ^"OO 



y 



dyy 2 

tc 2 Jo e v + 1 
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1(1 
2 16 



+ 1 

7T 2 7o e» + 1 ' 2vr 2 & 7o e« + 1 

1 f°° dyy 4 
T 3 27T 2 ^ /o e^ + 1 

_L_ Z^A. 

T 3^ + a T 3 e .' 



where 



T 3 r 00 dyy 2 
7r 2 7o e« + 1 : 



a = 



Jo 



ev+l 



Z J0 e a+l 



6.47, 



and thus, we get 



Pu 



rn ui h v + ^ 



T 



an,. 



i=l 



i=l 



2^ m^n„ II + I 



Pde = £V + V^(0 a ) = Y^rnuifiv U + 



+ V(<l>a) 



dp 



DE 



<9z 



E 

i=l 



oz 



dn u 
dz 



a \ A 2a <9£j 



+ 



^(00 

dz 



The stationary condition eq.(J2} leads to the relation: 

<9V(0 a ) ^ 9 



5 E 



5 E 



E ( X + ^) 



(41) 



(42) 



(43) 
(44) 
(45) 



i=l 



i=l 



12 



— n„ - 



dz d 
dz 



i=i 
dz 



5 E 



E 



<9mj, 7 - A a dn^ a 

-n,,-7T + m 



a 



i=i 



and thus, we have 
dz 



3 

E 

2=1 



dm L 



-n,, — a 



dz "e 



+ m i- 



9n w 1 



i ^6 



Using eqs.p3|l and (jUJ), the redshift dependence of the dark energy is 



#3* 



^DE 



(9,2 



" ' dm vi A <9n 



E 

i=l \ 
3 

+E 



i=l L 



dz 
dm. 



(72 J 

dm 



1 + _j_g m ^__ 



i=i 



m, 



dz 

dh u 
dz 



n,. — a 



A 1 dn^ 1 

dz Ci dz Ct 



9 fi 1^ 

£m V jTi v „ r. 



The equation of state parameter u> is 



w + 1 



[4 - 



3p 



DE 



where 



ft(T) 



i=i 



96 



E^(0 



Using eq. MT 



3 

i=l 



_ n„ 1 



E(^ + 



o 



3 a 

E & ~~ E 7~ 

i=i j=i si 

fc=l ;=i sZ 



n v 1 
T3^ 
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vi=l «, 



3 

fc=i 



i - 



E& 

i=i 



u'=i ?J , 



Thus, we have 



w + 1 



3 + 



E& 

i=i 



vi=i ^ 



2_^m vk n v + 2_^an u — I 

«Jfe=l i=i si/ 



3pde 

3 j. 2 ; :; 



3 E m ^™^ + 3 E + ~3 

i=l Si 



E f E 



i=l 



E& 



vj = l / fc=l 



3p 



DE 



2 / 3 a \ 3 T 



E& 

i=l 



E f E °^-f 

=1 S3 / k=l S* 



3p 



(E 



DE 



where the last term is negligible small because £i is much larger than 1. Thus we get 

U . T 



3 a x 3 



3 3 rj-i 

3 E m vi^v + 3 E afl vJ + ~3 ( E F I E m vkh v 

i=l i=l Si ™ \i=i sj / fc=i 



Pde(w + 1) 



3 

E m ^ 

1=1 



m ui n u + 

i=i \ 



3Ci 



Differentiating eq. (|53|) by the redshift parameter, we have 



<9ppE w dp BE 

dz dz 



" ; + FF" + ^e— = E 

i=l 



dw 1 

dz 



— — n u + m ui — 1 I 1 = I 

dz 



dz 3 \ dz & dz & £f 9z J 



(E 
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which leads to 



dp-DE dw 
— k — w + Pde^- 

dz dz 



3 

E 

i=l 
3 

E 

i=l 



-^n v + m vi — + ^[— lI +T --^ f - 



dp 



DE 



Then, using the relation (jlbj) . we have 



<9p 



DE 



(9,2 



( ' 8=1 



<9z 1/1 cfe 3 I £j £f (92 



9n w hah v /4T T % 



<9p; 



DE 



Oz 



dm,, 



+ 



dz dz ' 3 y & £f <9z 

Since the speed of sound squared in the dark energy is given by 

2 d P 

°s = a 

OpBE 

d(wp DE ) 

9pBE 
^DE„. ; , n ^dw 

9 PPE 

dz 

using eqs. (}4*%]) and (JSEJ), we get finally 



E 

8=1 



m h 



dn u 
dz 



E 

i=l 



dm 



dn u hah v ( 4T T <9£; 



<9z 



-n^ + m vl 



dz 



ii Q dz 



E 

i=i 



m,, 



dn u 

dz 



A dn v 

i=l 



dz 



Jt, dm ui 



5 A ' /5T Tgrry 



E 



m,. 



E 



<9,2 



<9z 
(55) 



(56) 



(57) 



(5f 



i=l ^ i=l 

It is easy to see that if the following relation is satisfied, the speed of sound squared becomes 
positive: 
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